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There is no tree-level flavor changing neutral current (FCNC) in the standard 
model (SM) which contains only one Higgs doublet. If more Higgs doublets are 
introduced for various reasons, the tree level FCNC would be inevitable except extra 
symmetry was imposed. Therefore FCNC processes are the excellent probes for the 
physics beyond the SM (BSM). In this paper, we studied the lepton flavor violated 
(LFV) decay processes h —>> (jlt and r /ij induced by Higgs -/i-r vertex. For 
r /i 7 , its branching ratio is also related to the htt, hr + r~ and hW + W~ vertices. 
We categorized the BSM into two scenarios for the Higgs coupling strengths near 
or away from SM. For the latter scenario, we took the spontaneously broken two 
Higgs doublet model (Lee model) as an example. We considered the constraints 
by recent data from LHC and B factories, and found that the measurements gave 
weak constraints. At LHC Run II, h /it will be confirmed or set stricter limit 
on its branching ratio. Accordingly, Br(r grf) < 0(1O -10 — 10 -8 ) for general 
chosen parameters. For the positive case, r —)> /iy can be discovered with 0(1O 10 ) 
r pair samples at SuperB factory, Super r-charm factory and new Z-factory. The 
future measurements for Br {h /it) and Br(r —>► fij) will be used to distinguish 
these two scenarios or set strict constraints on the correlations among different Higgs 


couplings, please see Table II in the text for details. 
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I. INTRODUCTION 


In the standard model (SM), we can diagonalize the gauge couplings and Yukawa cou¬ 
plings simultaneously, i.e., there is no flavor changing neutral current (FCNC) at the tree 
level. In the quark sector, flavor changing neutral currents occur at loop level with the help 
of CKM quark mixing matrix [lj. However, in the lepton sector, it is extremely suppressed 
by GIM mechanism j2] in the SM due to the smallness of neutrino mass. For example, for 
the lepton flavor violation (LFV) process t l —>■ % 7 

r(£ ->• tj'i) _ 3a 

r(£i —> ijVi&j) 327r 

in the SM |3j where V t] are the PMNS lepton mixing matrix [4j elements. With the data 
from neutrino oscillation [5j, it is estimated to be 

Br (fi ->■ e 7 ) ~ 0( HT 56 ) and Br(r -)■ e(p) 7 ) ~ e>(l(T 55 - 1(T 54 ) (2) 
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in the SM. It is far away from the recent experimental upper limit PC ! 1 

1.2 x 1(T 7 (Belle) f 4.5 x 1(T 8 (Belle) 

’ Br ( r < { 

3.3 x 10 8 (BaBar) I 4.4 x 10 8 (BaBar) 

and Br(/i —>■ e 7 ) < 5.7 x 10 -13 (MEG), all at 90% C.L. (3) 


Br(r —» e 7 ) < 


and the near future sensitivities with the improvement of an order MU. So that the 
discovery of the signals £ l —> (i.p at future colliders would clearly indicate new physics (NP) 
beyond the SM (BSM). Generally speaking the FCNC process will be one of the best probe 
of the BSM for future hadron and electron-positron colliders [12J. 

In July 2012, a new boson was discovered at LHC mm, and its properties are like 
those of a SM Higgs boson HE The Higgs mediated LFV process is attractive because of 
a 2.4 ct hint found by CMS Collaboration [1G| in the search for h —> /at process 2 . Assuming 
that the Higgs production cross section and total decay width are the same as those in the 
SM, the best fit (B.F.) branching ratio and 95% upper limit (U.L.) are respectively [16J 3 

Br (h ->• ht) = (0.84l°;p)% (B.F.) and Br (h pr) < 1.51% (U.L.) (4) 


1 For either B factory with L ~ 0.5ab 1 luminosity at ^/s = 10.6GeV (T(45) threshold). 

2 Recently the ATLAS Collaboration also published the searching result in the same process m with the 
result close to that in m by CMS Collaboration. 

3 For the full LHC Run I data with L cs 25fb -1 luminosity at ^s = (7 — 8)TeV. 
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If this signature was confirmed at future colliders, it would clearly indicate NP in the Higgs 
sector. In the extensions of SM, there may be direct Higgs- fx-r coupling to explain this 
hint, for example, in some types of two Higgs doublet models (2HDM) jTH] , like type III 
2HDM [T9H22] , 2HDMs with other flavor symmetries [231 - 125] , Lee model [26l 27], and other 
models [281430] . It may be also related to other phenomena like the excess in tth searches 
[31], b —> s semi-leptonic decays m, anomalous magnetic moment (g — 2) for // [32], LFV r 
decays [2D ESI E2HS3, or even the lepton flavored dark matter [HH]- Writing the Higgs-/i-r 
vertex as 


h 

£hfir ~\p2, 4"" Yth'TlIJ’R 3“ h.c.), 

and adopting the Cheng-Sher ansatz [36], the data gave [16] 


(5) 


^/|y /1T | 2 + |y T/1 | 2 <5xio- 3 or 


'(IVI 2 + |y T/1 |> 


2 m M m T 


< 2. 


( 6 ) 


In the future, at low energy e + e colliders like Super-B factory pun], Super r-charm 
factory [33138] or the new Z-factory [39], there would be signatures or stricter constraints 
for r —» /ry process; and at high energy colliders like LHC Run II at y/s = (13 — 14)TeV, 
there would be signatures or stricter constraints for h —>■ [at process. The results would 
be comparable and may give new constraints on the Higgs-yr-r coupling or the correlations 
among the couplings between Higgs and other particles. 


This paper is organized as follows. In |section H] we present the effective interactions 
and branching ratios for h —> /it and r —> py processes; |section ITT] and |section IV] contain 
the constraints from recent data and at future colliders respectively; |section V 
conclusions and discussions. 


are our 


II. EFFECTIVE HIGGS-p-r INTERACTION AND DECAY WIDTHS FOR 

h -> pr AND r —> py PROCESSES 


Based on 2HDM (type III), the higgs effective couplings can be written as 


= c v h [ ^W^W- + ^Z"Z a ) - h 


c t m t - c T m T _ 

tiJ-R H- t l t r + h.c. 


V 


V 


h 

71 


{Y^ulTr + Y t ^t l h r + h.c.) 


(7) 
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FIG. 1: Typical Feynman diagrams which contribute to r —/iy decay. The left one is Higgs- 
mediated one-loop diagram, the middle and right ones are W boson and top quark mediated 
Barr-Zee Type two-loop diagrams respectively. 

7 7 




where q stands for the coupling strength ratio compared with that in SM 4 and stands 
for the LFV coupling just like that in Eq. © With a direct calculation [16J, for the h —» /jt 
process, we have 

Br(ft^ M r) = I £^-(|y„ T | 2 + |r T „| 2 ) (8) 

where r/, means the total decay width of Higgs boson and in SM we have I\sm = 4.1MeV 
PTi] for rrih = 125GeV. 

The t —>■ fi'y decay process is loop induced. The dominant contribution usually comes 
from Barr-Zee type 5 two-loop diagrams since there is an additional ( m T /mh ) 2 log (rn\/m%) ~ 

Following the formulae in 1121. 


0(10 3 ) suppression in the one-loop amplitude [42], see the Feynman diagrams in Figure 1 


Br(r —» n'y) 
Br(r —¥ fivv) 


487r 3 « 

ri 2 

Lt p 


(l^l 2 + UkI 2 ) ■ 


(9) 


4 The c t and c T may be complex while cy must be real, and in the SM cy = c t = c T = 1. 

5 This type of two-loop diagrams was first proposed by Barr and Zee [41i during the calculation for lepton 
electric dipole moment (EDM). 
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Here the left (right) handed amplitudes Al(r) can be expressed as |2I3J 1,221 021 US] 6 


M*r) 

= "Al, 1-loop -loop) d - Al, 2-loop ("df#,2-loop) 


16^2 


: 7I" 
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3ivm T v 
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Crln (^) “ ^ Rc ( r d) - ylm(c r ) 
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Re(c t )/ ( r 4)+ilm(c t )g 

\ m hJ W 


where all the functions /, g and h come from 2-loop integrations as 


m = 5 
^ 1 


f 1 1 — 2x(l — x) ( x(l — x ) 

/ dx—--r-In - 

/ o x(l — x) — z \ z 


/ dx— --- 

l Q x(l — x) — z 
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»i 


h(z) = ~2 
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x{l — X ) 
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x(l — x) — z \ x{l — x) — z 
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m 


w 


mi 


( 10 ) 


( 11 ) 

( 12 ) 

(13) 


The small contributions from heavy neutral higgses, charged higgs and Z-mediated loop are 
all ignored. Defining 

(14) 

1 fJLT 1 T/J, 

the equation (|9j) should be changed to 

Br(r —>■ fi'j) 487t 3 ck|^4| 2 


Br(r —>■ \lvv) 


Lt p 


(IVI 2 + IW 2 ) 


(15) 


Here Br(r —» guu) = 17.4% from PDG |5j. 

For both decay processes, the LFV parameter comes in the form \/\Y lir \ 2 + |y r/1 | 2 , thus 
we do not need to study the details about the chiral properties of the LFV coupling. Since 
both Br oc (|1%-| 2 + \Y Tfl | 2 ), the ratio Br(r —>■ [h —>■ gr) does not depend on 

Therefore in this paper we will focus on the correlations among the Higgs couplings. 


6 The results in these papers are different. We checked the calculations and got the result consistent with 
that in [32 by Omura et. al. 
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FIG. 2: Distribution for R = Br(r —» /i 7 )/Br(/i /jlt) in cy — \ct\ plane in unit of (r/ ljtot )/MeV) 
fixing c T = 1. We take at = (0,7r/10, 7t/6) from left to right. The green regions are for R < 10 -10 ; 
the yellow regions are for 10 -10 < R < 10 -9 ; the blue regions are for 10 -9 < R < 10 -8 ; the cyan 
regions are for 10 -8 < R < 3 x 10 -8 ; the orange regions are for 3 x 10 -8 < R < 6 x 10 -8 ; the red 
regions are for 6 x 10 ' “ 8 < R < 10 7 ; and the brown regions are for 10 7 < R < 1.5 x 10 7 . 



0.8 1.0 1.2 1.4 1.6 1 . 


III. CONSTRAINTS BY RECENT EXPERIMENTS 


In general cases, a t = arg( q) and a T = arg(c T ) may be nonzero. The replacement 


<?h 


Br (h iir) —>■ — Br (h —>■ fir) = (cos 2 a t + 2.31 sin 2 a t )Br(h —> fir) (16) 

&h, SM 

should also be taken into account in 0 where cfh stands for the Higgs production cross 
section 7 and ct^sm means that in SM. To consider the numerical constraints on the couplings 
in ([7]) , we should take some benchmark points. Our fitting results [26j preferred |c r | ~ 1 
for almost all chosen for other parameters, so in this paper we take |c r | = 1. The regions 
cy < 0.4, \c t \ < 0.5 and \c t \ >2 are excluded for most cases by our fitting results, so we 
never consider those regions in this paper. 


According to (10), R = Br(r —> —>■ /it) is sensitive to the interplay between 

cy and c t . The cancelation between W loop and t loop induced amplitudes would make R 
very small in some regions especially for a t ~ 0. In |Figure 2| and Figure 3] we show some 
R = Br(r —> p 7 )/Br(/i —> /rr) distribution in cy — \c t \ plane in unit of (r^tot/MeV) for 
some different a t . From the figures, we can also see the cancelation behavior clearly when 
a t is small. For larger a t , the imaginary parts of the amplitudes would give more important 
contributions, and the imaginary parts of one loop contribution would also become more 


7 Gluon fusion process is dominant in this case. 
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FIG. 3: Distribution for R = Br(r —> p 7 )/Br {h — * /ir) in cy — |cj| plane in unit of (r/ ljto t)/MeV) 
fixing c T = 1. We take at — (vr/4, 7 t/2, 27t/ 3) from left to right. The green regions are for R < 10 -7 ; 
the yellow regions are for 10 -7 < R < 2 x 10 -7 ; the blue regions are for 2 x 10 -7 < R < 4 x 10 -7 ; 
and the cyan regions are for 4 x 10 -7 < R < xlO -6 . 



important as later figures shown. 

Here and in the following sections, we categorize BSM into two scenarios. In scenario I, 
we choose most Higgs couplings close to those in SM, especially cy ~ 1 and r^/r^sM- Since 
the experimental data [15] are consistent with the SM predictions, this scenario is popular. 
While the data still allow the Higgs couplings away from those in SM and these scenarios 
are attractive, because they are strongly related to BSM physics. In scenario II, we choose 
Lee model [261 EZj as such a benchmark model. Our previous work [26] showed that there is 
no SM limit for the lightest scalar in Lee model. We take the 125 GeV Higgs boson as the 
lightest one, so some of its couplings must be away form those in SM, especially cy should 
be small. In that paper, we considered full constraints by data and showed it is still alive. 
The fitting results for Higgs signal strengths allowed cy ~ 0.5, and at the same time, |q ,\ 
and F/, must be smaller than those in SM. The results are not sensitive to charged Higgs 
loop contribution. The typical T^/T^sm ~ (9(0.1) for different |c&| choice. In both scenarios, 
\c t \ ~ |c T | ~ 1 are preferred. 

In the scenario I, we take |c t | = 0.6,1.2,1.8 and plot the predicted branching ratios for 
t —> /i 7 in Figure 4 with cy = F/j/T^sm = 1 assuming Br {h —>■ /it) = 1.51% as the CMS 


upper limit, white regions are already excluded by recent data. For |c t | < 1.7, all the 
choices for ( a t ,a T ) are still allowed by recent data using this set of benchmark point, thus 
the recent r —> /ry measurements cannot give further constraints. While in the scenario 
II, the predicted branching ratios for r —» /ry are highly suppressed to be of 0( 10 -9 ) that 








FIG. 4: Br(r /iy) distributions in at — a T plane for for cy = T^/T^sm = 1, taking \c±\ — 
0.6,1.2,1.8 from left to right. The green regions are for Br(r /iy) < 1.5 x 10 -8 , the yellow 
regions are for 1.5 x 10 -8 < Br(r —>> fi 7 ) < 3.0 x 10 -8 , and the blue regions are for 3.0 x 10 -8 < 
Br(r fi 7 ) < 4.5 x 10 -8 . White regions are already excluded by recent data. 
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FIG. 5: Br(r —>► /iy) distributions in at — a T plane for cy — 0.5, r^/r^sM = 0.3, taking |<y| = 
0.6,1.2,1.8 from left to right. The green regions are for Br(r —>► fi 7) < 2.5 x 10 -9 , the yellow 
regions are for 2.5 x 10 -9 < Br(r /iy) < 5 x 10 -9 , the blue regions are for 5 x 10 -9 < Br(r —> 
/iy) < 7.5 x 10 -9 , and the cyan regions are for 7.5 x 10 -9 < Br(r —>► /iy) < 1 x 10 -8 . 





Lee model is not constrained by recent data. We take \c t \ =0.6,1.2,1.8 again and plot the 


predicted Br(r /iy) in Lee model in Figure 5 with cy = 0.5 and T^/r^sM — 0.3, assuming 
Br(/i —)> fir) = 1.51% as the CMS upper limit. 


IV. CONSTRAINTS AT FUTURE COLLIDERS 

Kopp and Nardecchia studied the phenomenology of h —» /it at future LHC (s/s = 
13TeV). With 300fb _1 luminosity, their results showed that for 07 , = cr^sM, if no signal is 
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observed, the expected upper limit at 95% C.L. should be set as Br (h —>■ fir) < 7.7 x 10 4 
which means 


Vl>VP + IV„P< 1.1x10- or 


< 0.45. 


(17) 


On the other hand, a signal would be observed at over 3 a if Br (h —>■ /it) > 1.3 x 10 3 which 
means 


+ !*,!»> 15X10- or 


> 0.59. 


(18) 


The SuperB factory is a e + e _ collider at T(4,S') threshold with the luminosity 75ab _1 . For 
the LFV decay r —> /ry, if no signal was observed at the SuperB factory, the expected upper 
limit at 90% C.L. should be set as (TOl 


(19) 


Br(r ->• /iy) < 2.4 x 10 9 . 


On the other hand, a signal would be observed at over 3cr if 


Br(r —> /iy) > 5.4 x 10 


-9 


( 20 ) 


At the Super r-charm factory, which is a e + e _ collider at y/s = (2 — 7)GeV with the 
luminosity 10ab _1 , there would be about 2.5 x 10 10 pairs of t + t~ [38] . And the sensitivity 
for LFV decay r —> /ry would be of (9(10 _1 °) [38] because of the suppression in background 
compared with that at SuperB factory 8 . And the same sensitivity (~ C7(10 -10 )) would be 
also achieved at new Z-factory [39J with 0(1O 12 ) Z bosons. 


For the r —>■ /ry results, there are three typical cases listed in |Tablc I| in which a positive 
result means a over 3a evidence and a negative result means an exclusion at 90% C.L. as 
usual. The typical choices are Br(r — > /iy) ~ 10 —8 ,10 -9 ,10 _1 ° for each case. And for h —> /j,t 


TABLE I: Choices for typical Br(r —>■ //,y) in different cases. 


The dominant backgrounds come from t + t ~7 events with a hard enough photon at SuperB factory; while 
at Super r-charm factory, \J7 is not far away above the t + t~ threshold that almost all photons from 
t + t _ 7 are soft. 
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Result at 

Result at 

Typical choice 


SuperB 

Super r-charm 

on Br(r —>► /iy) 

Case I 

Positive 

Positive 

~ 10“ 8 

Case II 

Negative 

Positive 

~ 10“ 9 

Case III 

Negative 

Negative 

< 2 x 10" 10 


results, we should consider the cases for LHC with positive or negative result separately. 


A. LHC with Positive Result 


A positive result in the h —» fir search would mean a direct evidence on LFV Higgs-/i-r 
coupling. We take (a/ l /cr/ li sM)Br(h —>■ /it) = 1.5 x 10 -3 , 3 x 10 -3 , and 6 x 10 -3 as benchmark 
points in this subsection. 

First, consider scenario I in section ITT] where the coupling strengths are close to those in 
the SM. Taking cy = = 1, |c*| = 1 and 1.5, we show the Br(r —» /ry) distributions 


in a t — Q't plane in |Figure 6] with the boundaries set according to the sensitivity of SuperB 


factory. We can see that if (cr/ l /a/ lj sM)Br(h —>■ fir) > (2 — 3) x 10 -3 , the typical predicted 
Br(r —>■ /ry) would reach the SuperB sensitivity. While if (a/ l /cr/ li sM)Br(h —> fir) was smaller, 
the t —y /xy process would not be found at SuperB factory. 

Then we should focus on the green regions which mean the cases with negative results 
at SuperB factory. Here we show the Br(r —> /ry) distributions in a t — a T plane in |Figure 7 
with the boundaries set according to the sensitivity of Super r-charm factory. For Br (t —> 
Pi) rsj 10 9 or smaller, \a t \ < 1.5 were favored. If LHC gave positive results, the typical 
predicted Br(r —>■ /r/) must reach the sensitivity of Super r-charm factory in this scenario. 
If Super r-charm factory gave negative results, it would give strict constraints on the Higgs 
couplings. 


In summary, For case I in |Table H if the SuperB factory gave positive results in searching 
t — y /jL'y (thus it must be discovered at Super r-charm factory as well), (a t ,a T ) would fall 


into the blue or cyan regions in Figure 6 The value of a T was usually free for larger | c t \ 
and (cpj/<Jh,sM)Br(h. —> fir), while a t \ > 1 were more favored for any case. While for case 
II in |Table I[ SuperB factory gave negative results but Super r-charm factory gave positive 
results, \a t \ < 1 would be favored, but for most cases there would be no constraints on a T . 
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FIG. 6: Br(r / 17 ) distributions in at — a r plane for for cy = r^/r^sM = 1, taking |q| = 1 
in the first line and |q| = 1.5 in the second line, and (ph/^h, sm)Bi*(/i —)• /it) = (1.5, 3, 6) x 10 -3 
from left to right. The green regions are for Br(r fi 7 ) < 2.4 x 10 -9 , the yellow regions are for 
2.4 xlCT 9 < Br(r —>• /J 7 ) < 5.4 x 10 9 , the blue regions are for 5.4 x 10 9 < Br(r —>• /ny) < 1 x 10 8 , 
and the cyan regions are for Br(r —7 // 7 ) > 1 x 10 -8 . 







For case III in Table I that both factories gave negative results, larger a T | and \c t \ would 
be favored. 


Second, consider Lee model which is scenario II in |section Hl| In this scenario, both cy 
and r fc /r ftiSM are smaller that the predicted Br(r —» fij) are smaller. For example, taking 
( c ’ 7 i./ cr /i,SM)Br(h -7 fir) = 3 x 10 -3 as a benchmark point, the predicted Br(r —» /i 7 ) < 
(0.8 — 1.6) x 10 -9 which cannot lead to a positive result at SuperB factory. 


We show the Br(r -7 /ry) distributions in a t — a T plane in Figure 8 with the boundaries 
set according to the sensitivity of Super r-charm factory, and all the colored regions are for 


Br(r —> n'y) < 2.4 x 10 -9 , thus case I in Table I would be disfavored. 

In this scenario, the results for Br(r -7- //yj cannot reach the sensitivity of SuperB factory 
but they will reach the sensitivity of Super r-charm factory. If Super r-charm factory gave 


negative results as case III in Table I, it would give strict constraints on the Higgs couplings 
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FIG. 7: Br(r / 17 ) distributions in at — a r plane for for cy = r^/r^sM = 1, taking |q| = 1 
in the first line and |<y| = 1.5 in the second line, and (<Jh/&h, SM)Br(/i —)• /it) = (1.5, 3, 6 ) x 10 -3 
from left to right. The green regions are for Br(r —)> /iy) < 2 x 10 —10 , the yellow regions are for 
2 x 10 -10 < Br(r —>► /iy) < 5 x 10 -10 , the blue regions are for 5 x 10 -10 < Br(r /iy) < 10 -9 , 
and the cyan regions are for 10 -9 < Br(r —>► /iy) < 2.4 x 10 -9 . 

Ctj- Ctj- QZj- 








as well that \a t \ < 1 would be favored but the constraints on a T would be weak. While if 


Super r-charm factory gave positive results as case II in Table I, larger \c t \ and a t would be 
favored. 


B. LHC with Negative Result 


In this subsection we choose (a/ l /cr/ lj sM)Br(h —>■ jir) = 7.7 x 10 -4 as the LHC expected 
95% C.L. upper limit together with the replacement (16). In scenario I in |section III| where 
the coupling strengths are close to those in SM, the predicted Br(r —>■ /17) < (1 — 2) x 10 -9 ; 
while in scenario II in section III, as the Lee model scenario, the predicted Br(r —y /ry) < 
(2 - 4) x 1(T 10 . 

We should discuss the two scenarios separately. We show the Br(r —> /./.y) distributions 
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FIG. 8 : Br(r —» /xy) distributions in at — ot T plane for for cy = 0.5 and T^/r^sM = 0.3, taking \ct\ = 
0.8 in the first line and |q| = 1.2 in the second line, and (<Jh/&h, SM)Br(/i —>► fir) = (1.5, 3, 6) x 10 -3 
from left to right. The green regions are for Br(r fij) < 2 x 10 — 10 , the yellow regions are for 
2 x 10 -10 < Br(r fij) < 5 x 10 -10 , the blue regions are for 5 x 10 -10 < Br(r fij) < 10 -9 , 
and the cyan regions are for 10 -9 < Br(r —/ 17 ) < 2.4 x 10 -9 . 





FIG. 9: Br(r —» /iy) distributions in ctt~a T plane for for cy = T^/r^sM — 1? taking |<y| = 1, 1 . 2 ,1.5 
from left to right. The green regions are for Br(r /iy) < 2 x 10 — 10 , the yellow regions are for 
2 x 10 -1 ° < Br (t —» /iy) < 5 x 10 _1 °, the blue regions are for 5 x 10 -1 ° < Br(r — » /iy) < 10 -9 , 
and the cyan regions are for 10 -9 < Br(r fij) < 2.4 x 10 -9 . 
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FIG. 10: Br(r —» /xy) distributions in cp — a T plane for for cy = 0.5 and T^/r^sM = 0.3, taking 
\ct\ = 1,1-2,1.5 from left to right. The green regions are for Br(r —> jin/) < 2 x 10 -10 , the yellow 
regions are for 2 x 10 -10 < Br(r —>• /xy) < 5 x 10 -10 . 





in a t — a T plane in |Figure 9] for scenario I and in |Figure 10| for scenario II respectively. If 
LHC gave negative results, the case I in |Table T| cannot appear thus we focus on case II and 
III. For scenario I, if Super r-charm factory gave negative results, a t | < (0.3 — 1) would be 
favored, else the other regions would be favored. For scenario II, most regions are allowed 


for case III in Table I that both e + e colliders gave negative results. 


V. CONCLUSIONS AND DISCUSSIONS 


In this paper we discussed the Higgs-/x-r coupling induced LFV decay processes h —» /xr 
and r —>■ /xy. For the later process, the branching ratio is also closely related to the htt, 
hr + r~ and hW + W~ couplings. We categorized the BSM into two scenarios, namely scenario 
I (II) with the Higgs coupling strengths close to (far away from) those in the SM, and for 
the latter scenario we took the Lee model as an example. We showed the possible numerical 


values of Br(r —»■ /xy) for different cases from Figure 4 to Figure 10 


If the future LHC run gives positive results on h —>■ /xr, different measurements on Br(r —> 
/xy) at super B factory and super r-charm factory would distinguish the two scenarios or 
imply the favored parameter choices. For case I in |Table I[ with positive results from both 
SuperB and Super r-charm factories, scenario I would be favored while scenario II would 
be disfavored or even excluded. For typical parameter choices, see the blue or cyan regions 
in |Figure 6] in details. For case II in |Table I[ with negative result from SuperB factory 
but positive result from Super r-charm factory, both scenarios would be allowed and some 
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constraints would be given on the Higgs couplings. See blue and cyan regions in |Figure 7 
and |Figure 8| for scenario I and II separately to find detail information on parameter choices. 
For scenario I, a T would be free for most cases, but regions near (a t , a T ) = (0, ±7r) would be 
disfavored for larger \c t \ and Br(/i —> /it) . For scenario II, | a t > 1 would be favored thus it 
implies large CP-violation in htt coupling. For case III in |Table~H with negative results from 
both SuperB and Super r-charm factories, scenario II would be more favored, but scenario 
I would not be excluded. See green regions in |Figure 7] and |Figure 8] for scenario I and II 
separately. 


If the future LHC run gives negative results on h —¥ /it, case I in Table I cannot be 


explained. If case I really happened, we would need other models. For case II in Table I 


scenario I with \a t \ > (0.5 — 1) would be favored, which implies large CP-violation in htt 
coupling. While there would be almost no constraints on a T . See |Figure 9] for details. For 
case III in |Table I[ nothing about LFV are to be seen at future colliders. Scenario I with 
l c t| ^ (0.5 — 1) would be excluded, while other regions for both scenarios are allowed. 

In |Table II| we summarize the implications corresponded to all the six future possibilities 
depending on the measurements of Br(/i —> /j,t ) at the LHC and Br(r —» ^7) at the super 
B factory and super tau charm factories. With the help of future measurements on LFV 
processes h — > /j,t and r — > /ry at both high and low energy colliders, for most cases, we 
would be able to distinguish different BSM scenarios or set constraints on Higgs couplings. 


P-IV in |Table II| would be strange. If it is really the case in the future, the Higgs induced 
LFV would not be the underlying reason. It would require other mechanism beyond Higgs 
sector to generate large enough LFV processes such as r —> /ry. 
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